МІНІСТЕРСТВО ОСВІТИ І НАУКИ УКРАЇНИ
НАЦІОНАЛЬНИЙ УНІВЕРСИТЕТ «ЛЬВІВСЬКА ПОЛІТЕХНІКА»
кафедра захисту інформації
Лабораторна робота №3
на тему: «ЗАСТОСУВАННЯ ХВИЛЬКОВОГО ПЕРЕТВОРЕННЯ
ДЛЯ ОБРОБКИ СИГНАЛІВ»
з курсу «Цифрова обробка сигналів та зображень»
Виконаd: ст. групи
ЗІД-12 Мазур Іван
Прийняв: Хома В.В.
Львів 2010
Мета роботи: ознайомитися зі змістом хвилькового перетворення, основними можливостями хвилькової обробки реалізованими в програмному пакеті MatLab, отримати навики розрахунку дискретного хвилькового перетворення, декомпозиції сигналу на апроксимуючі та деталізуючі складові, хвилькова очищення сигналів від шумів.
Завдання:
1. Ознайомитись з теоретичним матеріалом.
2. Скласти програму в середовищі MatLab згідно завдання.
2.1. Провести очищення сигналу від шуму з використанням хвилькового перетворення відповідно до завдання у табл. 2. Сигнал зберігається у файлі Lab_3_варіант у змінній signal.
Варіант
Тип вейвлету
Рівень
роз-кладу
Поріг
Тип
порогу
Масшта-бування
Сигнал
10.
Добеші 6-го порядку
4
sqtwolog
жорсткий
Тип 1
Lab_3_10.mat
Лістинг програми:
load Lab_3_1; % Завантаження сигналу
tn = (0 : length(signal)-1); % вісь часу
w_name = 'db6'; % тип вейвлету
Lev = 4; % рівень розкладу
type = ' sqtwolog '; % поріг
[Lo_D, Hi_D, Lo_R, Hi_R] = wfilters(w_name); % розрах-к фільтрів декомпозиції і реконструкції
% розрахунок АЧХ фільтрів
[h_D, f] = freqz(Lo_D, 1, 0:0.001:1, 2);
h1_D = freqz(Hi_D, 1, 0:0.001:1, 2);
[h_R, f] = freqz(Lo_R, 1, 0:0.001:1, 2);
h1_R = freqz(Hi_R, 1, 0:0.001:1, 2);
figure(1); plot(f, abs(h_D), f, abs(h1_D)); grid on; % графік АЧХ фільтрів декомпозиції
figure(2); plot(f, abs(h_R), f, abs(h1_R)); grid on; % графік АЧХ фільтрів реконструкції
[C, L] = wavedec(signal, Lev, w_name); % хвилькова розклад сигналу
cA4 = appcoef(C, L, w_name, Lev); % добування коефіцієнтів апроксимації сА3
[cD4 cD3 cD2 cD1] = detcoef(C, L, [4 3 2 1]); % добування коефіцієнтів деталізації cD4, сD3, cD2, cD1
% Побудова графіків коефіцієнтів розкладу
figure(3); subplot(511); plot(cA4); grid on; xlim([0 L(1)-1]); ylabel('cA4');
subplot(512); plot(cD4); grid on; xlim([0 L(2)-1]); ylabel('cD4');
subplot(513); plot(cD3); grid on; xlim([0 L(3)-1]); ylabel('cD3');
subplot(514); plot(cD2); grid on; xlim([0 L(4)-1]); ylabel('cD2');
subplot(515); plot(cD1); grid on; xlim([0 L(5)-1]); ylabel('cD1');
A4 = wrcoef('a', C, L, w_name, 3); % Обчислення апроксимуючої складової А4
D4 = wrcoef('d', C, L, w_name, 4); % Обчислення деталізуючої складової D4
D3 = wrcoef('d', C, L, w_name, 3); % Обчислення деталізуючої складової D3
D2 = wrcoef('d', C, L, w_name, 2); % Обчислення деталізуючої складової D2
D1 = wrcoef('d', C, L, w_name, 1); % Обчислення деталізуючої складової D1
L_sig = length(signal); % Розмірність сигналу
% Побудова графіків сигналу та його часових складових
figure(4); subplot(611); plot(signal); grid on; xlim([0 L_sig-1]); ylabel('S(n)');
subplot(612); plot(A4); grid on; xlim([0 L_sig-1]); ylabel('A4');
subplot(613); plot(D4); grid on; xlim([0 L_sig-1]); ylabel('D4');
subplot(614); plot(D3); grid on; xlim([0 L_sig-1]); ylabel('D3');
subplot(615); plot(D2); grid on; xlim([0 L_sig-1]); ylabel('D2');
subplot(616); plot(D1); grid on; xlim([0 L_sig-1]); ylabel('D1');
% Оцінка рівня шуму деталізуючих коефіцієнтів
for ii = 1 : Lev
stdc(ii) = wnoisest(C, L, ii);
end;
% Обчислення порогу для кожного рівня
for ii = 1 : Lev
thr(ii) = thselect(detcoef(C, L, ii)/stdc(ii), type);
thr (ii) = thr(ii) * stdc(ii);
end;
% жорстке порогування коефіцієнтів cD1
for ii = 1 : leng...